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Abstract 

Consider the following evolution model, proposed in |Hj by Bak and Sneppen. Put N vertices 
on a circle, spaced evenly. Each vertex represents a certain species. We associate with each 
vertex a random variable, representing the 'state' or 'fitness' of the species, with values in [0, 1]. 
The dynamics proceeds as follows. Every discrete time step, we choose the vertex with minimal 
fitness, and assign to this vertex, and to its two neighbours, three new independent fitnesses 
with a uniform distribution on [0, 1]. A conjecture of physicists, based on simulations, is that 
in the stationary regime, the one-dimensional marginal distributions of the fitnesses converges, 
when N — * oo, to a uniform distribution on (/, 1), for some threshold / < 1. 

In this paper we consider a discrete version of this model, proposed in In this discrete 
version, the fitness of a vertex can be either or 1. The system evolves according to the following 
rules. Each discrete time step, we choose an arbitrary vertex with fitness 0. If all the vertices 
have fitness 1, then we choose an arbitrary vertex with fitness 1. Then we update the fitnesses 
of this vertex and of its two neighbours by three new independent fitnesses, taking value with 
probability < q < 1, and 1 with probability p = 1 — q. We show that if q is close enough to 
one, then the mean average fitness in the stationary regime is bounded away from 1, uniformly 
in the number of vertices. This is a small step in the direction of the conjecture mentioned 
above, and also settles a conjecture mentioned in [2J. 

Our proof is based on a reduction to a continuous time particle system. 

Key words: species, fitness, evolution, interacting particle system, self-organised criticality, 
coupling, contact process, stationary distribution. 
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1 Introduction 

The Bak Sneppen model, introduced in [3], has received a lot of attention in the literature, see for 
instance [T], [S], [Ej, jl], [H], P and [12]. In [T], it is described how Bak and Sneppen were 
looking for a simple mathematical model which was supposed to exhibit evolutionary behaviour, and 
which was also supposed to fall into the class of processes showing self-organised critical behaviour. 
For physicists, self-organised critical behaviour refers to power law decay of temporal and spatial 
quantities, without fine-tuning of parameters. After many attempts, Bak and Sneppen arrived at 
the following process. 

Think of a system with N species. These species are represented by N vertices on a circle, evenly 
spaced. Now each of these species is assigned a so called 'fitness', and in this model, the fitness is a 
number between and 1. The higher the fitness, the better chance of surviving the species has. The 
dynamics of evolution is modelled as follows. Every discrete time step, we choose the vertex with 
minimal fitness, and we think of the corresponding species as disappearing completely. This species 
is then replaced by a new one, with a fresh and independent fitness, uniformly distributed on [0, 1]. 
So far, the dynamics does not have any interaction between the species, and does not result in an 
interesting process. Interaction is introduced by also replacing the two neighbours of the vertex with 
lowest fitness by new species with independent fitnesses. This interaction represents co-evolution of 
related species: if a certain species becomes extinct, this has an effect on other species as well. The 
neighbour interaction makes the model very interesting from a mathematical point of view. 

It is extremely simple to run this model on a computer. Simulations then suggest the following 
behaviour, for large iV (see [K] and P for simulation results). It appears that the one-dimensional 
marginals are uniform (in the limit for N — ► oo) on (/, 1) for some / whose numerical value is 
supposed to be close to 2/3. This threshold value / is the basis for self-organised critical behaviour, 
according to [3], PP and jS]. Since in the limit there is no mass below /, one can look at so called 
avalanches of fitnesses below this threshold: start counting at the moment that there is one fitness 
below / and wait until all fitnesses are above / again. The random number of updates, for instance, 
counted this way, is suppose to follow a power law, and there is no fine-tuning of parameters. 
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It is a challenge to prove any of the above statements. Note that in order to prove power law 
behaviour, one should first prove the existence of the threshold / with the property that in the 
limit for iV — > oo, all one- dimensional marginals are concentrated on (/, 1). Indeed, one can define 
avalanches corresponding to other thresholds as well, but it is not expected that these avalanches 
have power law behaviour. This is only expected (and observed) for the self-organised threshold /. 
Therefore, this should be the starting point of a rigorous mathematical analysis of the model. 

Simple as the model may appear, it turns out to be very difficult to say anything at all about 
the limiting one-dimensional distributions. It is therefore natural to try to prove a similar result in 
a simpler model. In this light, we have chosen to study a discrete version of the model, which was 
proposed in 2J, and which can be described as follows. Fitnesses of species can now only be or 1. 
The dynamics in this simpler model proceeds as follows: at every discrete time step, we choose an 
arbitrary vertex with fitness 0. If there is no such vertex, then we choose an arbitrary vertex with 
fitness 1. We update the fitnesses of this vertex, and of its two neighbours, by three new independent 
fitnesses, taking value with probability < q < 1, and 1 with probability p — l — q. This process is 
called the BS process in this paper. We show that if q is close enough to one, then the mean average 
fitness in the stationary regime is bounded away from 1, uniformly in the number of vertices. This is 
a small step in the direction of the conjecture mentioned above, and answers a question which was 
posed in j2|. 

It should be noted that this discrete version of the model does not show self-organised critical 
behaviour. Nevertheless, we think that understanding of the discrete model also increases our under- 
standing of the original model, if not in a technical sense, certainly in a conceptual sense. Admittedly, 
the discrete version suggested here and in |2j is only one out of many possible discrete versions, but 
we see no reason to complicate matters unnecessarily by choosing more complicated discrete versions. 
The reader will notice that the proof of our main result is already quite complicated. 

In order to state our main result, here follows some notation. As before, the number of vertices 
is denoted by N, and we denote by r)isr(n)i the state of the z-th vertex after n updates of the process. 
We will prove the following result. 
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Theorem 1.1 If q is close enough to one, then there exists c q > such that for any N G N and 

ie{l,...,N}, 

lim P{r, N {n) t = 0) > c q . (1.1) 

Note that lim P{rj^{n)i = 0) exists, because 77jv(n) is a finite state, irreducible and aperiodic Markov 
chain. 

In the next section, we reduce the problem to a problem in a continuous time, monotone particle 
system. In this system we will be able to prove results uniformly in N, by exhibiting graphical 
representations and an infinite space version of the particle system. 



2 Reduction to a monotone continuous time process. 

In this section we define a useful continuous time stochastic process £(£), independent of N. We 
construct the process £(t) via a graphical representation. The graphical representation GR is a 
random graph on the space-time diagram Z x M + . We define GR via a set of independent so called 
bundles {II fc } fceZ , where each bundle n fe consists of eight independent Poisson processes on R, 

n, = {nr,nr,...,nr,n^}, 

with parameters g 3 /(l — g 3 ), q 2 p/(l — q 3 ), ■ ■ ■ ,p 2 q/(l — q 3 ),p 3 /(l — q 3 ) respectively. (We use the factor 
1/(1 — q 3 ) to rescale time in a convenient way, as will become clear later.) For each process n^ 1 ' " 2 ' " 3 
we perform the following procedure. At i-th arrival t j ^' <72 ' <T3 of Yl CT j }' a2 '' 7 ' i , j 6 Z, we draw arrows in 
ZxR from {k,r^' a2 ^) to (k - 1, r^'" 2 '" 3 ), iff a x = 0, and from (k^f 2 '* 3 ) to (k + 1, r^'" 2 '" 3 ), iff 
o- 3 = 0. We draw a * in Z x R at every (k+j, r^' a2 ' a3 ) with ^ = 1, j = -1, 0, 1. We say that (x, tj) 
is connected to (x, £2) by a time segment, if ti < £2 and there are no *'s on (x, t), £ G [ti, ^2)- An open 
path'f^)] 1 ! 2 with trace (xo, so), • • • , (x n , s n ) is a map from [£i, £2] to Z such that t\ = s < ■ ■ ■ < s n = £ 2 , 
j(t) = Xi, for t G [sj,Sj + i), < i < n, 7^2) = ^n, and every pair (xj,Sj), (xj + i,Sj+i) is connected 
either by a time segment or by an arrow. For any finite B C Z, x,y £ B, and t\ < t 2 G K, write 
(x, £i)~^(y, £2) in if there exists an open path 7(£)|i 2 in GR with 7(£i) = x, 7^2) = y, and 

7(t) G £, for all £ G [t 1? t 2 ]. 
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Figure 1: The graphical representation GR. (0,0) is connected to (— by an open path. The trace is (0,0), 
(0,r„7), (~l,r 7), (-l,T°i° ), (-2,T°\° ), (-2,r°_%), (-l,r°°° ) (-l,t). 

For any finite iCBcZ and i, s > we denote by C^f'^s) the random subset of vertices x G B 
such that there exists y = y(x) G A, and (y,t)-^(x,t + s) in GR\b- 

The definition of £s ( s ) via the existence of certain open paths implies a number of useful 
properties. The first is monotonicity: 

^*\s) C s > if i C C, B C D, (2.1) 

The second is the semigroup property: 

(a) $^(*) = {0}, (22) 

(b) zt BA ' t){si) ' t+Sl \s 2 -s l ) = & A > t \s 2 ), HACB,0< Sl <s 2 . 

The monotonicity property (|2.1j) allows us to define the process £g *^(s) for any i C B C Z: 

= Km dr'^). 

B'Tfl, (2.3) 

s' finite 

Note that due to the monotonicity property (j2.1j) . the limit at the r.h.s. is independent of the sequence 
B' | Z. The process is now defined as follows: 

^) = ef 0) W, f>0. (2.4) 
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We now extract the BS process 7/7v( n ) from the graphical representation GR as to have r)N{n) 
and defined on the same probability space. 

The N vertices are labeled by A(N) = {-N' -1, . . . , N" + 1}, where N' + 3 + N" = N, N" — 1 < 
N' < N", and the observation site is labeled by 0. We define and r(i) to be the left and right 
neighbours of i, respectively, with appropriate boundary conditions: 



l(i) 



i-1, ie[-N',N" + l], 
N" + 1, i = -N' - 1, 



ru 



T "\, 



i + l, i G [-N' - 1, N" 
-N' - 1, z = iV" + 1. 

A state of the BS process is determined by the subset of the vertices in state 0. Thus the state 
space Sn of the BS process consists of the all subsets of A(N). If we denote the state of a site i in a 
configuration r\ G by rji, then we have an identity 



Vi 



l{i£v}- (2-5) 



This identity is natural because the 0's play the 'active' role in the dynamics of BS process. It is 
possibly also slighty inconvenient for mathematicians, who are used to work with subsets of sites in 
state 1, (like in the contact process or in the 1-dim sendpile model, ect.). Indeed, for those processes 
we would have rji — l|i G T]^. Nevertheless, we will not reverse the roles of l's and 0's, and will 
work with ([2.5)1 . because of the conventional definition of the BS-model. 

We will now extract from GR two independent sequences of random variables U = (Uj), V = (Vj), 
and then we will define the BS process in terms of those sequences. Let II(AQ denote the superposition 
of all the Poisson processes associated to the vertices in A(A r ), i.e., with abuse of notation, 



n(A0= |J (n™unru.--uni 1() un 



in 

k 



fceA(V) 

Then II(AQ is a Poisson process on R with intensity N/(l — q 3 ). Denote by 



r(JV) = (r 1) r a ,...) (2.6) 
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the arrivals of IT(iV) after time zero. For every j 6 N there exists, with probability one, a unique 
Uj G A(iV) and Vj G {0, l} 3 such that Tj is the arrival of 11^.. It is clear that U, V and t(N) 
are independent and each consists of i.i.d. random variables. Note that Uj, j G N is uniformly 
distributed, that is, P(Uj = i) = 1/N, i G A(iV). The sequence U will be used as a sequence of 
random vertices-canditates for the update procedure. The distribution of Vj is simply the joint 
distribution of three independent Bernoulli random variables, taking value with probability q and 
1 with probability 1 — q. The sequence V will be used to determine the states of the vertices after 
the updates. 

We will define ^(n) inductively via the (random) increasing sequence (j n ) C N. Let jo = 0, 
Vn(0) — 0- Let n G N, and suppose that j n -i and r] N (n — 1) are already defined. If r] N (n — 1) = 0, 
then j n = j n _i + N + 1, rj N {n) = {Uj n }, i.e. we skip N elements of the sequences U and V, and then 
restart our process from the site Uj n . The reason to skip iV elements is that we want to have the 
following property: the more particles in state 1 we have in r]N(n — 1), the more elements of U and 
V, in mean, we skip to define ^(n). We will use this property later, in the proof of Lemma [2. II If 
VN{n — 1) 7^ 0, we wait until we choose a vertex in state 0: 



j n = min jj > j n _i | Uj G rj N (n - 1) j. 



VN(n)i 



and change the state of site Uj n and its neighbours l(Uj n ) and r(Uj n ) according the value of Vj n 
(0-1,0-2,0-3), i.e. 

' <Tl, 1 = l\Ujn), 

02, i = U jn , 

0-3, i = r(u jn ), 

rj N {n — l)j, otherwise. 
This finishes the construction of the process r)N(n). 

We will now introduce an 'intermediate' continuous time process £,§(t): 



r] N {0), tG[o,n), 

VN(n(j)), t G [Tj,T j+1 ),j > 1, 
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where (n(j)) is defined as 

n(j) = max jn £ N : j n < j j. 

It is clear that £,§(t) is a continuous time Markov chain on Sn- Observe that the processes £jy and 
?7tv are related via a random time change. If there are many l's around, then we typically skip more 
steps, so l's tend to be preserved in £jy. This intuition is articulated in the following lemma. 



Lemma 2.1 We have 



lim P(v N (n) = l)< lim PU§(t) = 1 
PROOF: We prove ()2.7|) in two steps: 

(a) lim P[r] N (n) = 1 ) < lim P[r) N (n(j)) = 1), 

(b) lim P(^(n(j)) = l) = lim PU§(t) = 1 



We prove (b) first. We write 

oo 

P(f£(t)o = 1) = ^P(^(n(j)) = 1, and r,- < t < r j+1 ). 



j=0 



The random variable i] N (n(j)) is independent of Tj and Tj +1 . Hence 



^(eJv(*)o = 1) - lim P(vN(n(j)) = 1 
£ { P (w(n(j'))o = l)P(r, < t < r j+1 )} - lim p{ji N (n(J))o = 1 

oo 

£ ^(^("0*))o = l) - lim P(^(n(j))o = l) 



x P{ T j <t < T j+i) ~ * 0) as t oo, 



becauce r,- — ► oo, as j — ► oo in probability. This proves (b). 
For fa), we write 



(2.7) 



(2i 



p fc = lim P y~] r] N (n)i = k , 

n^oo \ ' I 

. ieA(JV) 
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and 



It is then clear that 



\ieA(N) 



P(VN(n) = 1) = £ *** 



and 



AT 

%W)). = l)=Ej' 

fc=0 

Now observe that when there are k < N vertices with fitness 1, the number of trials before we select 
a vertex with fitness has a geometric distribution with parameter (N — k)/k, hence the expected 
number of trials is equal to N/ ( N — k) . It follows that for k < N, 

N 



N-k 



Pk 



and similarly 



qN = 



(N + l)p N 



EJIo 1 W=iPi + (N + l)p 



N 



We write a k = k/N, b k = N/(N - k) for k < N, and b N = N + 1. Then 



P(VN(n(j)) = 1) = Y,~W 

k=0 



E^o h iPi 

Efc=0 a kbkPk 
E^O 6 lPl 

Now observe that for any probability vector p , . . . ,p^, and non-decreasing sequences < a < ■ ■ ■ < 
ajv and < b < ■ ■ ■ < b^, we have 

N \ / N \ N 

^a-kPk ) I ^2b k p k J < ^2a k b k p k , 

,fc=0 / \fe=0 / fc=0 

which can be proved by induction. Applying this general fact to the a k s, b k s and p k s above, we find 
that the last quotient is bounded below by J2 k=0 a k p k which is just P(r] N (n) = 1). (Note that here 
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we have used the fact that we skip iV choices if all vertices have fitness 1: this makes the sequence 
(bk) increasing.) This proves (a). 

□ 

The next step in the proof of Theorem 11.11 is the following lemma, where we relate the process 
£jv(£) to the graphical representation GR. To simplify notations further we will write GRn instead 
of GR \ [_jv',iv"] ■ 

Lemma 2.2 For any t x ,t 2 > 0, x,y G [-N',N"\, if £§(ti) x = and (x,t 1 )**+(y,t 2 ) in GR N , then 

t§(h)y = 0. 

Proof: Let (x, ti)^(y, t 2 ) in GR N , with trace (x , s ),..., (x n , s n ), i.e. every pair (x u Sj), (x i+1 , s i+ i) 
is connected either by a time segment or by an arrow. The statement will follow by induction, if 
we prove that Xi G £jv( s i) implies x i+ i G £^(s i+ i), for every < i < n — 1. Let x« G £,§{si). 
If (xj,Sj) is connected to (xi + i,s i+ i) by a time segment, then there are no symbols V on (xi,t), 
t G [si,s i+ i). Hence, there are no arrivals at the time interval [sj,s i+ i) at n^ ,1,<73 , n^ 1 ^ 2 ' 1 and 
n^,^' 73 , (o"i, cr 2 , C3) G {0, l} 3 . Thus, = Xi G £jj(s i+ i), because only the above arrivals can delete 
Xi from £jy(si+i). If Xi + i = Xi + 1 and (xj, s,) is connected to (xj+i, Sj+i) by an arrow, then there is an 
arrival at n£ 1,<72 '°, for some 01, 02 G {0, 1} at time Sj = Sj+i, and hence, G £§(si+i). Similary, if 
x i+ i = Xi — 1 and (x^, Sj) is connected to (x^+i, by an arrow, then there is an arrival at n°'. <J2,<73 , 
for some a 2 , cr 3 G {0, 1} at time = s i+1 , and again, x^+i G 

□ 

The last two lemmas imply that we have reduced the problem to GR. Therefore, in the next 
section, we work in this graphical representation. 

3 Proof of Theorem 11.11 

To simplify our notation £b ( s )> we will skip the upper index (A, t), if A = B and t = 0. We will also 
skip the lower index B, if B = Z. So for example, we will write ^ A,t ^(s) instead of £2 ^(s), Cb,oo)(^) 
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instead of Cr^^) \t), and C(-oo,x](*) instead of The ^ ea ^° cou Pl e by a graphical 

representation the processes with various lower indices is taken from [T3j . 

Note that for any t > 0, with probability one, £[_jv',oo)(£) 7^ and £(_oo,iv"] (*) 7^ 0- Thus we can 
define l-N'(t) and rjv»(t) as the leftmost and rightmost 0's of the processes £[-n> ,oo)(£) and £(-oo,iV"](X)j 
respectively. The following lemma was inspired by inequality (5.2) in 

Lemma 3.1 

(a) Z[- N >,oo)(t) 2 [l-N>(t), oo) n £(*), t > 0, 

(b) £(-oc,N»](t) d (-oo,nv»(*)] n^(t), i > 0. 

PROOF: Let y G [Z_jv/(t), oo) fl £(t). Then there exists X\ G Z and an open path 7i(s)|q in GR such 
that (sex, 0) is connected to (y, t) by 7i(s)|g. Since l-N'(t) G £[-jy,oo)(£) there exists 22 G Zfl [— iV', oo) 
and an open path 72(5) |q, laying completely within [— N', 00) x R such that (22, 0) is connected to 
(L N ,(t),t) by 

72 ( s ) I o- Let s * ^ e defined as 

s* = inf{s>0 : 72 (s) < 7l (s)}. 

Note that by definition any open path is a cadlag function of time, thus 72(5*) < 71 (s*). If s* = 
then the open path 7i(s)|q lays completely within [— N', 00). If s* > 0, we define the open path 

7s(s)lo as 

, , J 72OO, se[o,4 

[ 7i(s), s G [a*,*], 

see Figure 121 The open path 73(5) | has endpoint (y,t) and lays completely within [— N',00) x R. 
Thus G £[_jv,oo)(*)- 

The statement in (b) can be proved in a similar way. 

□ 

Define [I, r] N (t) as 

[l,r) N (t) =£[-N',oo)(t) D £(_«,,#'/](*), t > 0. 
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Figure 2: The open path connecting (x 2 ,0) to ( 72 (s*),s*) to (y,t) lays completely within [-iV',oo) x R. 

Lemma 3.2 If [l,r] N (t) ^ 0, for all t G [£i,t 2 ], then for any y G [^rjjvt^) there exists an open path 
7 (i)|g with 7 (t 2 ) = y and 7 (*) G [Z,r]*(t), t G [Ma]- 

PROOF: Let A' be the subset of t* G [ii,t 2 ] which have that for all y G [l,r]iy(t*), there exists an 
open path 7 (s)|** with 7 (t*) = y and 7 (t) G M^O), t G [Ml- If [Z,r] w (t) ^ 0, for t G [Ma] then 
l- N '(t),r N (t) G [-iV',iV"], for all t G [Ma]- Thus the event t* G A' is determined only by the GR 
inside [-N',N"\ x [Ma]- 

Suppose t* G A' and r is the first arrival of n(iV) after time t*. Then due to the definition of an 
open path [t*,r) G A'. We will prove that r G A'. Then, by induction, t% G A', because there are 
only finitely many arrivals in II(iV) at time interval [t\, £2] and t\ G A'. Let y G [I, t}n(t). Then there 
exist open paths 7 i(t)|t 1 G [— N', oo) and T2 C^) i ti e ( — °o,iV"] with endpoint (y,r). Since l~N'(t*) 
is the leftmost point of £[-jv,oo)(**), we have 7i(t*) > l^ N >(t*). Similarly, 72 (i*) < r N "(t*). The 
distance | 7 i(t*) — 7 2(^*)| < L because there is only one arrival of IT(iV) in the time interval (t*,r]. 
It follows from the above that 7 i(£*) G [l,r]jsr(t*) or 7 2(t*) G [Z,r]jv(t*). Suppose 7 i(t*) G [£,r]jv(£*)- 
Then 7i (t) G [l,r] N (t), t G [t*,r]. Since t* G A' there exists an open path 7 (s)|£ with 7 (f) = 7 i(t*) 
and 7 (t) G [i, r]jv(t), i G [ii,i*]. Hence, the path 7 3(i)|£" , defined as 



73 




7 (t), *e[M*), 
7l (t), te[t*,r}. 



has endpoint (y,r), and satisfies 73 (i) G [Z,r]jv(^), £ G [ii,T]. See Figure|3]for an illustration. Thus, 
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Figure 3: At the realisation [/,r]jv(r) = [l,r] N (t*) = {y}. Since (y,t*) is connected to (y,r) by 7(t)|[. and 4* e A', 
there exists an open path with endpoint (y, t), laying completely within [Z,r]jv(t), £ € [ti , 7-] . 



We need a number of results which are very similar to the corresponding results for the standard 
contact process. Therefore, we omit the proofs of the following three lemmas. Their proofs are 
modifications of the proofs of Theorem 3.19, Theorem 3.21 and Corollary 3.22 in ^U|. Note that in 
these three lemmas, we require q to be close enough to one. It is at this point where the rescaling of 
time by a factor 1/(1 — q 3 ) comes in. With this rescaling of time, the intensity of arrows gets large 
when q gets close to one. The (omitted) proofs of the three lemmas involve comparison with oriented 
percolation, and in order to make sure that the appropriate oriented percolation model percolates, 
we need a high intensity of arrows. 

Lemma 3.3 If q is close enough to one, then there exists u(q) > such that for any t > 



□ 




(3.1) 
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Lemma 3.4 Let 16Z and let l x (t) be the leftmost zero of £r X) oo)(£)> an d r x (t) be the rightmost zero 
of £(-00,2] If <? is close enough to one, then there exist C\(q), C2(q) > 0, depending only on q, such 
that for any m G N and t > 5m 2 



Pyx(s) > x + m, for some s E [t — 5m 2 , t] J < Ci(g)e ca fa') m ) 

/ \ (3-2) 

-Pf^z(s) < a; — m, for some s £ [t — 5m 2 , t] J < c\(q)e~ C2 ^ m . 

Lemma 3.5 For q close enough to one there exists c(q) > 0, depending only on q such that 

P (3f G [0, JV/2] : (0,0)->(JV,O in Gi? [0jOo) ) > c(g). 

Now we partition the time axis into intervals of length N, and we call the interval [iN, (i + l)iV) the 
z-th level. 

Definition 3.6 We call level i normalising if there exist t,t',t" G [iN, (i + 1)N) and x G [—N',N"] 
such that £§(t) x = 0, (x,t)^(x — N,t') in G-R(_oo,x] and (x,t)^(x + N,t") in GR[ X>00 ). 

See Figure 0] for an illustration of this definition. 



(i+l)N 



(x-N,f) 





N 



(x,t) 



iN 



-N' 



N 



Figure 4: The i-th level is normalising. 



We will use normalising levels to connect different open paths in GR. For this to work, we have to 
make sure that there are enough normalising levels. This is the content of the following key lemma. 
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Lemma 3.7 For q close enough to one, N large enough, and T > N 2 + N, there exists c(q) > 0, 
depending only on q, such that the probability to find no normalising level among levels \T/N\ — 
N,..., [T/N\ is at most e~ c ^ N . 

Proof: The events 

ji-th level is normalising j, i G N 

are not independent. But we will construct independent events which guarantee that certain levels 
are normalising. To do this carefully, let Ti(N) be the cr-algebra generated by the restriction of GR 
to Z x [0, iN + N). Let 

t* = min ji/V + N/2, inf{* > iN : £%(t) H [-JV', iV"] ^ 0}}, 

and 




min{x G [-JV 7 , iV"] : £§(t*) x = 0}, if t* < iN + N/2, 

0, if t* = iN + N/2. 



Let Ai(N) be the event that (x*, t*)-^(x* + N, t') in GR [xhoo) , and (x?-iV,t") in G J R ( _ 00jX * ] , 

for some G + iV/2). Note that Aj(A^) n {t* < iN + N/2} implies that the z-th level is 

normalising. Since t* is a stopping time and because of Lemma 1331 there exists p q< A > 0, depending 
only on q, such that 

P(A{N)) >p q , A , (3.3) 

uniformly in N. Observe that Ai(N) is ^-i(iV)-measurable, and is independent of J-i~\(N). 

We are now going to make sure that {t* < iN + N/2} occurs often enough. Let s* be defined as 
the smallest element s of the i-th level for which one of the following three conditions is satisfied: 

(1) $00 = 0; 

(2) £§(s)n[-N',N"]^<b; 

(3) s = iN + N/4. 

Note that s* is a stopping time (with respect to the natural filtration) and that s* < iN + N/4. We 
will now give a condition in terms of GR within the i-th level which ensures that s* < iN + N/4. 
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If £jv(iiV) = or £§(iN) n [— iV', iV"] 7^ 0, then s* = iiV. The remaining cases are those where 
is not empty and a subset of {— N' — 1,N" + 1}. We now pretend that £§(iN) is not empty 
and a subset of {— N' — 1,N" + 1}, giving three possible situations, namely £§(iN) = {—N' — 1}, 
£§(iN) = {N" + 1} or £§(iN) = {—N' - 1, N" + 1}. In each of these situations, we compute 

inf js > : £%{iN + s) = or £*(iN + s) n [-N', N"] ^ 0}, 

which we denote by S 1 , S 2 and S 3 respectively. We denote the maximum of these three numbers by 

S*: 

S* = max{S\S 2 ,S 3 }, 

and define the event B^N) as 

Bi(N) = {S* < N/4}. 

Note that Bi(N) is measurable with respect to the a-algebra generated by GR restricted to the 
i-th level, and therefore the Bj(iV)'s are mutually independent for different i's. Also observe that 
occurrence of B^N) implies that s* < N/4, and that 



lim Pi B 



i(N)) = 1. (3.4) 



N^oo 

Next, define Ci(N) as the event that the (N + l)-th arrival of H(N) after time s* takes place before 
time s* + N/4. Since T1(N) has intensity of order N, the number of arrivals of n(AT) in a time interval 
of length Af/ 4 has a Poisson distribution with mean of order A^ 2 . This implies that 



lim P(Ci(N) = 1. (3.5) 

Finally, D^N) is defined as the event that the first arrival of T1(N) after time s* is at a vertex in 

[-N', N"}, or that there is no such arrival during the time interval [s*, s* + N/4). It is clear that 

lim p(Di(N)) = 1. (3.6) 

The events Ci(N) and D^N) are A") -measurable, and independent of jF i _ 1 (A r ). Now we have 
that 

{^(AO n a(A0 n A(JV)} c {t* < *iv + iv/2}, 
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and therefore 

^Ai(N) n Bi(N) H Ci(N) n A(JV)) C ji-th level is normalising J. 
Also, for N large enough, we have, according to (|3.3j> . (|3.4jl . ()3.5|) and ()3.6|) that 

p(A(iV) n Bi(iV) n Ci(JV) n A(JV)) > ci(g), 

for some Ci(g) > 0, uniformly in iV. We may now write, using the independence of all events 
P( none of the levels \T/N\ — N, . . . , [T/iVj are normalising 



H /iV j 

p( p| (A-WnB^jnCi^nAW)' 



LT/7VJ 

< 

i=LT/JVJ-JV 
LT/7VJ 

II / , ((A(iV)ns i (iV)nC i (JV)nA(JV)) £ 

i=[T/AfJ-Af 



for some c(g) > 0. 



□ 



Proof of Theorem ll.lt Due to the symmetry of BS-process we can work with i = 0. According 
to Lemma 12.11 it suffices to prove that for q close enough to one, there exists c q > such that for 
any iV sufficiently large, 

lim P(£#(t) = 0) > c q . 

For any T > and iV we have 

p(t§(T) = o) > P(e(T) = 0) - p($(T) = 1, e(T) = o). (3.7) 

The first term in ()3.7j) is independent of ./V and positive, according to Lemma f3. 31 Thus it is remains 
to prove that 

p(^§(T) = 1, f (T) = o) -> 0, as JV -> oo, uniformly in T > T(A r ), (3.8) 
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for some T(N) < oo. Let T > N 2 + N. According to Lemma f3 .41 and the stationarity of GR, for q 
close enough to one, there exists c\ = ci(q), C2 = 02(g) > depending only on q, such that 



(i) P(l- N >(t) > 0, for some t E [T — N 2 — N, T] \ < c x e 

(ii) P(r N „(t) < 0, for some t G [T - N 2 - N, T]j < c x e~ 

It follows from ()3.9|) and Lemma l3~7l that 



-c 2 N 
c 2 N 



/ 



< p 



p(^(T)o = 1,£(T) = 



^(T) = 1, f (T)o = 0, 



Z_jv(t) < < r^(t), for all £ G [T — N 2 — N,T], 
y 3i G {|T/NJ - N, . . . , [T/N\} : i-th level is normalising y 



+cie 



-c 2 iV 



(3.9) 



(3.10) 




Figure 5: The main idea of the proof. With probability close to one, {£(T) = 0} implies {£jv(T) = 0}. 
Consider the first term of ()3.10|) . It follows from Lemma f3. II that the event 



l-N'(T) < < r N »{T), £(T) = 
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implies that {0} G [1,v]n(T). Then, by Lemma f3. 21 the event 

{0}e[l,r] N (T), 1 
l~N'{t) < r N „{t), for all t G [T - N 2 - N, T), J 

implies that there exists an open path 7(i)|y_ JV 2_ Ar with endpoint (0, T) and laying completely within 
[—N',N"\. If the 2-th level is normalising then, by definition, there exists t,t',t" G [iN, (i + 1)N), 
x G [-N',N"] and open paths 7z(s)|*' : {x,t)«*(x - N,t') in GR^^, j r (s)\f : (x,t)«*[x + N,t") 
in GR[ Xt00 ). Either 7/(s)|f or 7 r (s)|*" intersects r ){t)\lj,_ N 2_ N -, see Figure Hence (x, t) is connected 
to (0, T) by an open path in GR^, and by Lemma l2~2l we have £^(T) = 0. Hence the first term at 
the right hand side of (j3.1(J|) equals zero and the second term gives us the theorem. 

□ 

Acknowledgement: We would like to thank Jeff Steif and Jeremy Barbay for helpful discussions. 
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